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We have been looking into the vibrational distribution 

of O 2 in the atmospheres of Venus and Mars to compare them 

. + 

to that found on the Earth. Vibrationally excited 0 2 on 
Venus is formed in charge transfer reactions and lost in 
dissociative recombination. The vibrational levels are in- 
terchanged in quenching reactions and fluorescent scattering. 
In order to compute the vibrational distribution, we needed 
to compute the excitation rates and transition probabilities 

for fluorescent scattering in the second negative system of 
x 9 2 

°2 n g^ ' Because the potential curves for the X and 

-f 

A states of 0,, have minima at very different internuclear 

4 . 

distances, excitations from the X state can populate very 

high vibrational levels of the A state, which may then decay 

to high vibrational levels of the X state. Published values 

of the transition probabilities included only about ten 

levels of each state. This was completely inadequate for 

our purposes, so we recomputed the wavefunction and energies 

2 

for 61 vibrational levels of the X state, and 34 vibra- 

2 

tional levels of the A state, and the oscillator strengths 
for all of the transitions between these levels. We then 
computed the transition probabilities for the A-+X transitions, 
and excitation rates for absorption of a solar photon using 
our oscillator strengths and the solar fluxes of Vernazza 
et al. (1976) . {See attached preprint) . 
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We have calculated the vibrational distribution of 
°2 Ng) and °2 ^ n u^ ; *' n t * ie da y time Venunian ionosphere 
for models which are based on Pioneer Venus data. We find 
that at low altitudes quenching limits the fraction vibra- 
tionally excited to less than 5%. At high altitudes the 
vibrational distribution approaches that assumed for the 
reaction: 

0 + + co 2 0 2 + (v) + CO . (1) 

We have assumed that the reaction proceeds by the formation 
of a collision complex so that the available energy is 
statistically distributed among the degrees of freedom of 
the complex. This is probable for low energy encounters. 

Thus we have assumed that all of the available vibrational 
levels are populated equally in (1) and the relative popula- 
tion of each of vibrational levels 0 to 5 of 0 2 approaches 
.16 at high altitudes. 

The fraction of C> 2 vibrationally excited at the exobase, 
about 200 km, depends on the assumed value for the rate coef- 

* 4 " ■ 

ficient for quenching of 0 2 (v) by atomic oxygen: 

0 2 + (v) + 0 -*■ 0 2 + (v-1) + 0 (2) 

The rate coefficient for reaction (2) has not been measured. 

We first assume that the reaction proceeds with a rate coef- 
ficient of lxl0~ 10 cm 3 s"'' 1 ' and that the vibrational quantum 
number changes by only one in each reaction. For this model, 
49% are vibrationally excited at the exobase. This will 
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affect the exothermicity of the reaction: 

C> 2 + + e 0 + 0 

both because the average vibrational energy is near .' l .4eV 
at the exobase and because the dissociative recombination 
may proceed through different channels for different vibra- 
tional energy levels {Guberman private communication, 1983). 
Thus the hot oxygen corona will be affected by the vibrational 
distribution of 0 2 + (v) at the exobase and therefore the hot 
hydrogen corona will also be affected by scattering with 0 
atoms. (See attached preprint of a paper given at the XXVth 
COSPAR meeting and submitted for publication in Advances in 
Space Research) . 
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0 2 + is the major ion in the lower ionosphere of Venus, 

\ 

even though 0 2 is not an important neutral constituent (Taylor 

•f . 

et al* 1980). 0 2 is formed in the reactions 

1 4 - 

0 + C0 2 4 CO + °2 (v) * 1 * 33eV * d) 

and x 

0 + + C0 2 4 CO + ° 2 + ( v) + 1 , 20eV . (2) 

Reaction (1) is the major source below and reaction (2) above 
150km. The ions produced in these reactions will probably be 
vibrationally excited. Reaction (1) can populate up to v=6. 

For reaction (2), vibrational levels Jap to v=5 are energetically 
allowed. 0^ is lost mainly by dissociative recombination 

0 2 + <v) + e 0( 02 3 P, 1 D) + 0( 3 P, 1 D ; 1 S) (3) 

* 

The potential .-curves published by Guberman (1983) show that 
the repulsive state of 0 2 through which the dissociative recom- 
bination occurs will vary with the vibrational level of the 
0 2 + . Thus the rate of reaction (3) , the electronic states of 
the 0 atoms produced and hence the total exothermicity will 
also depend on the vibrational level of the ion. Guberman 
(private communication) has found that the O(^S) yields for 
0 + in the v=l and v=2 states are considerably enhanced over 
that for v=0. If the production rates of 0( 1 S) and 0( 1 D) are 
affected, • the airglow intensities of lines originating in 
these levels will also be affected. 
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Abreu et al. (1983) have examined data from the visible 
airglow experiment on Atmosphere Explorer-E and have found 
that the quantum yield of OC^'S) (0.09-0.23) varies with the 
ratio of the electron density to the atomic oxygen density. 

+ 

This ratio is a measure of the vibrational excitation of 0 2 
in the terrestrial ionosphere (Bates and Zipf , 1981) . 

Dissociative recombination of 0 2 + is the major source 
of hot oxygen atoms in the upper atmosphere of Venus (Nagy et 
al., 1981). The energy distribution of these atoms will depend 
on the exothermicity of reaction (3). One mechanism for the 
escape of hydrogen from Venus is collision of H atoms with 
hot 0 atoms (McElroy et al., 1982; Kumar et al., 1983). The 
predicted escape rate will therefore depend on the channel 
assumed for reaction (3). 

In this paper, we calculate the vibrational distribution 
of 0 2 + in the ionosphere of Venus for a model based on data 
from the Pioneer Venus neutral mass spectrometer (Hedin et al. , 
1983; Neimann at al., 1980). This model is essentially the 
same as that described by Fox (1982a, b) . The calculation of 
the ion densities includes both chemistry and diffusion. The 
relative populations of the vibrational levels are computed 
assuming photochemical equilibrium and are valid to about 230km. 
We included vibrational levels v=0 to v=60 of the ground state 
X 2 n g of 0 2 + , and v=0 to v=33 of the A 2 IT u state. 
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The X 2 n^ State 

+ . 

The sources and sinks of 0 2 in the Venusian ionosphere are 
summarized in Table 1. Reactions (1) and (2) far exceed photo- 
ionization and electron impact ionization of 0 2 as sources of 
4 “ + 2 

C >2 . 0 { D) may also react with C0 2 , although only the rate 

+ 2 

coefficient for charge transfer of fast 0 ( D) has been measured 

+ 2 

(Moran and Wilcox, 1978) . 0 ( D) has been found to react at 

near gas-kinetic^ rates in collisions with several molecules 

(cf. Glosick et al., 1978; Rowe et al. , 1980). We assume that 

the reaction occurs with a rate coefficient of 8x10 ^cm^s \ 

near that for N 2 and CO, and that the products are CO + 0 2 , 

+ 4 

as they are for reaction of 0 ( S) . The vibrational distribution 
+ 

of the 0 2 produced in these reactions is unknown. We assume 

that the available levels are all populated with equal probability. 

As sources of 0 2 (v=l) , reactions (1) and (2) are of 

similar magnitude up to 150km. Above that altitude, reaction 

(2) dominates. Dissociative recombination is the major loss 
+ 

process for C> 2 over the entire altitude range of interest. 

The value of the dissociative recombination coefficient will 
vary with the vibrational level of the ion, but no measurements 
or calculations of the magnitude of the. difference exist. We 
therefore assume that the rate coefficient measured by Mul and 
McGowan (1979) applies to all the vibrational levels. 

In addition to the processes which cause net production 
or loss of 0 2 + , there are processes which interchange the vibra- 
tional levels. They are presented in Table 2. Fluorescent 
scattering of sunlight in the second negative system of C> 2 
is such a process. We have computed new excitation rates and 
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transition probabilities from the dipole moment computed by 
Wetmore et al. (1984). For solar fluxes, we used the values 
presented by Vernazza et al. (1976). 

Vibrationally excited 0 2 + (v>0) may be quenched in colli- 
sions with atoms or molecules. The rate coefficient for 
quenching by C0 2 and several other molecules has been reported by 
Bohringer et al. (1983) for v=l and v=2. We assume that the 
rates for higher levels are equal to that for v=2. 0 2 (v>0) 
may also be quenched by atomic oxygen (reaction 12) . Our stand- 
ard model assumes that the rate coefficient is lxlO~^^cm^s ^ 
for all v, but we also explore the consequences of larger or 
smaller values. In these reactions, we assume that the vibra- 
tional level of the 0 2 changes by only one quantum, as one 
would expect for vibration-translation (V-T) energy transfer 
(Bohringer et al., 1983) 

At low altitudes, quenching is more important as a source 

or sink for v~l than the reactions which cause net production 
+ ... 

or loss of C> 2 . The vibrational distributions at low altitudes 
will not directly reflect the distributions produced in reactions 
(1) or (2). 

The computed relative populations of vibrational levels 
v=0-2 and v> 3 for altitudes from, 115 to 250km are shown in Fig. 1. 
Below 140km quenching limits the fraction vibrationally excited 
to less that 5t>. At high altitudes, the populations of vibrational 
levels 0 to 5 approach 0.16 , since the major C> 2 production 
mechanism, reaction (2), is assumed to produce each of those 
levels with equal probability. The fraction found in the ground 
state decreases gradually above 140km. The rate of this reduction 


5 


+ 

depends on t..e assumed value for quenching of 0 2 (v) by atomic 
oxygen (reaction 12). Figure 2 also shows altitude profiles of 
the fraction found in v=0 for quenching coefficient k^ 2 of 
lxl0~^cm 3 s"”^ and 6x10 ^°cm^s A quenching coefficient of 

GxlC'^cm^s** 1 should be considered an upper limit, since this 
value is essentially gas kinetic. The results of Abreu et al. 

(1983), however, probably exclude values this large. 

Table 4 gives the fractional population in v=0 to 9 at 205km 
for each of the models considered. In the standard model, 59% 

are vibrationally excited at that altitude. If k^ 2 is reduced to 

-11 3-1 . . , 

1x10 cm s , the fraction vibrationally excited increases to 

79%; if k^ is increased to 6x10 ^cm^s ^ fraction in v>0 decreases 
to 22%. We also consider a model in which the vibrational level 
produced in reactions (1) and (2) is always the highest level 
possible. This assumption should give an upper limit to the frac- 
tion vibrationally excited. Table 5 shows that for this assumption, 
90% of the 0 ^ are in v>0, and for v=l to 5, the fraction in each 
level increases with the vibrational quantum number. 

* « » 4 * 

Fluorescent scattering in the second negative system of 0 2 

is important as a source of v>6. The only other production mecha- 

, . + 2 

nism for these levels is reaction of 0 ( D) with C0 2 , which may 

produce 0 2 (X H^) in vibrational levels up to v=24 and possibly 

+ 4 4 

0 2 (a n ) up to v-3, although production of the a Jl u state would 

involve a spin change. The second negative band system is weak 
and spread'over a large wavelength range. Figure 2 shows the in- 
tegrated overhead intensities summed in 50A bins from 1800A to 
900 0A. The total intensity is about 640R. 
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Discussion 

The vibrational distribution of 0 0 on Venus is significantly 

M 

different from that on Earth (Fox and Dalgarno, 1984). On Earth, 
the major production mechanisms are photoionization and electron 
impact ionization and the reaction 

0 + + 0 2 + 0 2 + (v<7) +0 . 

The vibrational distribution of 0 2 + produced in ionization of 
0 2 is determined by Franck-Coudon factors and is different from 
the statistical distribution assumed for the charge transfer 
reactions. The presence of 0 2 as a quencher is also important 
for Earth. In general quenching of a molecular ion by its parent 
neutral is efficient. Because the reaction occurs by symmetric 
charge transfer, the change in v is not limited to one quantum 
number. The final state of the 0 2 + (v) will be influenced by 

-f. 

Franck-Condon factors for 0 2 0 2 . At high altitudes on Earth 

+ ... 
the distribution of C> 2 (v) approaches 50% vibrational ly excited. 

Both the rate of dissociative recombination of 0 2 (v) and the 

final electronic states of the 0 atoms produced will be influenced 

by the vibrational level of the ion. The major source of both 

O^S) and 0(^D) in the Venusian thermosphere^ is photodissociation 

of C0 2 (Fox and Dalgarno, 1981} . The peak of this source occurs 

at altitudes (115km) at which 0(^D) atoms are quenched. The total 

dissociative recombination rate is comparable to or exceeds the 

production of 0(^D) due to all other sources for altitudes above 

140km. Thus we expect the branching ratio for 0 ( D) production 
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* 4 * , 

in dissociative recombination of 0 2 to influence the intensity 
and lineshapes of the 6300A 0(^D -*■ 3 P) emission. 

0( S) is not as readily quenched and th© density of this 
species peaks near 115km. The dissociative recombination rate 
exceeds the production of 0(^S) only for altitudes above 130km. 
Thus, although we do not expect the totaJ overhead intensity 
of the emissions from O(^S) at 5511k end 2972A to be greatly 
affected by the value of the yield of O(^S) in (3), the volume 
emission rate above 130km will reflect the yield if it is signi- 
ficant. 

Nagy et al. (1981) predicted the presence of a hot oxygen 
corona around Venus analogous to the terrestrial oxygen corona 
(Yee et al. , 1980). The major source of hot oxygen atoms was 
found to be dissociative recombination of 0 2 . The axothermici ty 
of the reaction and therefore the energies of the 0 atoms produced 
will depend on the vibrational state of the 0 2 + ion. At 205km, 
the average vibrational energy is ,38eV in the standard model. 

"3 1 

If the product atoms are in the 0( P) and 0( D) states, 
the energy released is 5.0eV, so the addition of 0.38eV will not 
be significant. If, however, the products are 0( 1 S)+ Ol^D), 
the exothermicity is only 0.83eV for C> 2 + (v=0). The addition of 
one or two quanta of vibrational energy may significantly change 
the energy distribution of the 0 atoms produced. Killeen and 
Hays (1983) examined line profile measurements of the oxygen 
green line in the terrestrial ionosphere and concluded that the 
0( 1 S) + 0( 1 D) channel is favored over. the 0( 1 S) + 0( 3 P) channel 
by a factor of four. This is in agreement with the theoretical 


t 



3 

prediction of Guberman (1983) that the major channel for 

production of 0( 1 S) is through the repulsive state of 0 2 

1 1 

which dissociates into one 0( S) and one 0( D) . If so, the 

line profile of the 5577A emission should reflect the vibra- 

• + 
tional state of the 0 2 ion. 

More impprtant to the total exothermicity however, is 

the possibility .-that the reaction proceeds through different 

channels for different initial vibrational levels of C> 2 . 

McElroy et al. (1982) proposed that the escape of hydrogen 

from Venus is caused by collisions of thermal H atoms with hot 

0 atoms produced in 

0 2 + + e + 0( 1 D) + 0( 3 P) + 5. leV . 

Kumar et al. (1983) have evaluated the source and found that 
it is important in the current atmosphere, but that charge 

4» t 

exchange of hot H with H dominated in the past. Table 3 shows 
that near the exobase, the 0 2 + is considerably vibrationally 
excited. If recombination in different vibrational levels 
proceeds through different chanrals, the energies of the oxygen 
atoms will change and therefore the hydrogen escape rate due 
to this source will be affected. 
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Figure Captions 


1. Relative populations of 0 2 in v=0, 1, 
over the altitude range 1 X S to 250km. 

2. Integrated overhead intensities of the 

negative bands of 0 2 + (A 2 n u -X 2 n_) . The 
are summed in 50A bins. •' 
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Table l. Sources and Sinks of in the 

Venusian thermosphere 




Reaction 


Rato (cm 2 s”^) 

Reference 

1 

0 + C0 2 + -+■ 0 2 + (v<6) 

+ 

CO 1.64xl0~ 10 

Fehsenfeld et al. 
1970 

2 

0 + + C0 2 -► 0 2 + {v<5) 


CO 9.4xlQ“ 10 

Smith et al . 1978 
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1978 
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Smith et al. 1978 

6 
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Fehsenfeld 1977 
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Goldan et al, 1966 
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c -*• C0‘ +0 
-* C* + 0, ■ 


5X10’ 11 

5xl0 -11 

Prasad and Huntress 
1980 


Table 2. Processes that interchange the vibrational levels of 
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Table 3, Vibrational distributions of at 205km in the 

Venusian ionosphere for several models. The numbers 
given are the fraction in each vibrational level. 

Model : 

vibrational 1234 5 

level 


0 

.41 

.20 

.78 

.10 

.61 

1 

• 15 

.17 

.064 

.059 ■ 

.062 

2 

.14 ^ 

.17 

.054 

.090 

.068 

3 

.12 

.16 

.043 

.14 

.075 

4 

.098 

.15 

.030 

.22 

.082 

5 

.062 

.13 

.017 

.35 

.090 

6 

.0053 

. 0057 

.0019 

.017 

.0048 

7 

.0033 

.0022 

.0014 

.0050 

.0017 

6 

.0027 

.0017 

.0012 

.0047 

.0015 

9 

.0021 

.ooi: 

.0010 

.0036 

.0012 

10 

~.01 






Model: 

1. Standard model: k 12 = IxlO -10 ernes’* 1 . 

-11 3 -1 

2. reduced to 1x10 cm s 

3. k^ 2 increased to 6x10 10 cm^s 

4. Production of O- in reactions (1) and (2) is 
assumed to populate the highest vibrational levels 
possible. 

. Same as 4 but with k^ 2 increased to 6xl0 - ^ cm 3 s - ^. 
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